Chromosome Research

Recent technological advances allow for the large-scale acquisition of genetic and morphological data: high-throughput sequencing has transformed the field of genomics while synchrotron X-ray microtomography enables rapid, noninvasive 3D imaging. However, integrating these approaches for the same specimens is challenging because X-rays can fragment DNA, and DNA extraction damages internal morphology, particularly relevant for small bodied organisms, such as insects. We systematically tested multiple extraction protocols and irradiation conditions across three model insect species. We irradiated more than 1,000 specimens under varying conditions and tested DNA quality through DNA barcoding and UCE sequencing. Our results demonstrate that high-quality DNA and high-resolution tomograms can be obtained from the same individuals, provided that the parameters are carefully optimized and rapid SR-CT scanning precedes DNA extraction. In this respect, our findings establish practical guidelines for combining genomics and phenomics, paving the way for comprehensive integrative digitization of biodiversity.

Whole genome duplications leave lasting traces in our genomes. How these present in terms of gene content and order varies over time. While collinear blocks of paralogs, long stretches of conserved gene order and content termed microsynteny, are a distinctive feature of comparatively recent WGD and have been integral in reconstructing the history of ancestral duplication events, this signal degrades over time, making analysis of older events non-trivial. While gene order degrades quickly, gene content is often better conserved and recent work takes advantage of this to reconstruct older events and ancestral pre-WGD and post-WGD chromosomes. However, these new methods are complicated and not well-documented. Here we develop an automated and user-friendly pipeline for reconstructing ancestral chromosomes before and after WGD, and use the conservation of gene content to infer chromosomal rearrangement events in this timeframe. We verify the efficacy of our tool by reconstructing the ancestral acipenseriform, a model system for vertebrate WGD and rediploidisation. Our pipeline should serve to make ancestral reconstruction more accessible and provide a solid foundation for future analysis.

Whole-genome duplication (WGD) causes chromosome instability through multipolar chromosome segregation driven by supernumerary centrosomes. WGD cells formed through distinct processes, mitotic slippage (MS) and cytokinesis failure (CF), show a prominent difference in viability after multipolar chromosome segregation: MS causes a more skewed homologous chromosome distribution than CF, resulting in more frequent nullisomic chromosome segregation with poorer survival through the first mitosis. However, the determinants of route-dependent differences in post-WGD cell viability remain largely unknown, particularly regarding the contribution of spatial rearrangement of supernumerary centrosomes. Here, we found marked differences in supernumerary centrosome distribution upon entry into the first mitosis after MS and CF, stemming from distinct nuclear geometry. The distinct centrosome distributions differentiated kinetochore capture patterning after MS and CF, whereas their modulations had minimal effect on the fidelity of subsequent chromosome segregation. In contrast, artificially extending the centrosome-centromere capture range by depleting the microtubule depolymerizer MCAK drastically suppressed the MS-linked aggravation of nullisomic chromosome segregation through equalizing chromosome capture by each supernumerary centrosome. These results suggest that centrosome-centromere capture range, rather than the spatial arrangement of the centrosomes themselves, determines the fidelity of chromosome segregation after WGD. Our findings provide fundamental insights into atypical cell proliferation mechanisms after WGD.

Dinoflagellate chromosomes adopt a highly condensed and organized morphology, with periodic bands and arches observed by traditional Transmission Electron Microscopy (TEM). However, the limited two-dimensional field of view of TEM has prevented a precise characterization of the inherently three-dimensional organization of dinoflagellate chromosomes. Moreover, given the vast diversity among dinoflagellate species and the lack of a systematic comparison of their chromosomes, it remains unclear whether dinoflagellate chromosomes share common organizational features or instead exhibit significant cell- or species-specific differences. Here, we acquire three whole-nucleus 3D Focused Ion Beam Scanning Electron Microscopy (FIB-SEM) datasets at 4 nm voxel size for each of four dinoflagellate species: Symbiodinium microadriaticum, Breviolum minutum, Fugacium kawagutii, and Crypthecodinium cohnii. We compile these data with previously published image volumes from four additional species and present an analysis of the largest collection of dinoflagellate FIB-SEM images to date. Common features observed across all eight species include the absence of physical confinement or spatial clustering of chromosomes in the nucleus. In addition, by decomposing each chromosome into a weighted sum of orthogonal shapes using Spherical Harmonics Expansion, we find a principal component encapsulating 88% of the total shape variance that is common to all species. However, our analysis also reveals differences in chromosome morphology across species. First, while many chromosomes exhibit surface ridges with left-handed helical twist, the proportion of chromosomes with such ridges varies extensively across species. Second, while chromosomes in most species are discrete and well-separated, chromosomes in F. kawagutii are interconnected in a single contiguous network. Lastly, to our knowledge, we report the first observation in eukaryotic cells of toroid-shaped DNA objects, whose numbers vary dramatically across cells and species. Overall, our results show that dinoflagellate chromosomes exhibit both shared organizational features and pronounced species-specific deviations.

Satellite DNA is long arrays of tandem repetitive DNA located often near the centromeres of chromosomes, whose function, or lack of, has been debated since its discovery. Although situated in heterochromatin, satellite DNA may be expressed as long noncoding RNAs (lncRNAs). Although there are a few examples of satellite lncRNAs being characterized, and functions suggested, how widespread and functionally important they may be for developmental processes is not understood. Here, we take an evolutionary approach to investigate satellite lncRNA expression in Drosophila spp. ovaries, a tissue whose development is well-characterized but where satellite expression has only been minimally explored. Using a publicly-available total RNAseq dataset, we find that 118/156 surveyed satellite DNAs were expressed across 10 species, with 33 satellites having high expression over 20 RPM. However, all but two of these expressed satellites (AAACTAC in D. virilis and ACAGACAGACAGG in D. ananassae) had higher read counts in a sister smallRNA dataset, suggesting that most satellite transcripts primarily serve as precursors for piRNA biogenesis. The two "stand-alone" lncRNAs were highly strand-biased, with 96-97% of the total reads coming from one strand. We further investigated AAACTAC expression with RNA FISH and found the transcript is specifically present in the oocyte nucleus following a dynamic spatiotemporal pattern, with the highest expression in stage 3-5 oocytes. The transcription pattern of AAACTAC is conserved in the three other virilis clade species that contain this satellite DNA. Further, we found expression of unrelated satellites in more distantly related D. borealis and littoralis both in the oocyte and the nurse cells. Overall, our work identifies a novel lncRNA AAACUAC found in the early oocyte nucleus, which is conserved across ~5 MY of evolution, and is therefore a strong candidate for the discovery of novel functions of satellite lncRNAs in development.

Bipolar spindle assembly is essential for accurate chromosome segregation. KIFC1, a conserved Ran- regulated minus-end-directed kinesin-14 motor, accumulates in the nucleus during interphase and promotes chromatin-mediated spindle assembly during mitosis and meiosis. In human oocytes, reduced KIFC1 levels destabilize meiotic spindles, a defect that can be rescued by increasing KIFC1 expression. However, how KIFC1 expression levels affect mitotic spindle stability during cleavage divisions in vertebrates remains unclear. Here, we show that whereas an approximately 50% reduction in KIFC1 causes no detectable defects in spindle assembly, approximately 10-fold overexpression of KIFC1 induces monopolar spindle formation, leading to chromosome mis-segregation and embryonic lethality in medaka early embryos. KIFC1 overexpression results in ectopic centrosomal localization during interphase, impairing the separation of duplicated centrosomes before mitotic entry. Analyses of KIFC1 mutants demonstrated that these centrosome separation defects require KIFC1s microtubule-binding and motor activities and are further enhanced by deletion of KIFC1s nuclear localization sequences. Together, our findings demonstrate that tight regulation of KIFC1 expression and its nuclear sequestration is essential for the proper separation and positioning of duplicated centrosomes before mitotic entry, thereby ensuring efficient bipolar spindle assembly during the rapid cleavage divisions of vertebrate embryos. HighlightsO_LIKIFC1 accumulates in the nucleus and at the embryonic spindle midplane via the Ran pathway. C_LIO_LIPartial KIFC1 depletion does not impair spindle assembly in medaka early embryos. C_LIO_LIKIFC1 overexpression induces monopolar spindles by preventing centrosome separation. C_LIO_LICentrosome separation defects require KIFC1 microtubule-binding and motor activity. C_LI

Plumage pigmentation is under intense sexual and natural selection and plays an important role in the speciation process in birds, so there is much interest in uncovering the genomic basis of plumage colour differences between populations and species. Three species of North American woodpeckers, the red-breasted (Sphyrapicus ruber), red-naped (S. nuchalis), and yellow-bellied sapsuckers (S. varius), provide a particularly promising opportunity to unravel the genomic mechanisms of plumage colour differentiation. The breeding ranges of the three species are mostly non-overlapping but adjacent, with hybrid zones occurring where the ranges meet. The species pair with the most similar plumage colouration (S. varius and S. nuchalis) is not the most closely related pair genomically (S. nuchalis and S. ruber is), providing an opportunity to determine the subset of the genome that underlies the plumage colour variation. Using admixture mapping of whole genome sequences from each species and hybrids from each species pair, we show close associations between the colour of multiple plumage patches and wide swathes of the Z-chromosome. These results highlight how comparable changes in one sex chromosome can cause either slight plumage pigmentation changes (S. varius vs. S. nuchalis) or large-scale shifts from dimorphism to monomorphism and from primarily black and white to primarily red plumage colouration (S. varius and S. nuchalis vs. S. ruber).

Cellular systems must act robustly to maintain organismal health, including maintaining biophysical properties that allow for appropriate cellular function, and adapting these properties through changes such as those that occur during cell division. However, we still lack tools to measure many of these physical properties with precision in the living cell. For example, the mechanical properties of the nucleoplasm, the fluid-like substance that fills the nucleus, have not been fully characterized. To investigate these properties, we have turned to the fission yeast Schizosaccharomyces pombe (S. pombe), a well-established, genetically tractable model organism that has been used extensively for studying a variety of cell biophysical processes and structures, including the cytoskeleton and cell division. It is an apt system for studying how the nucleus adapts over the course of the cell cycle, since it undergoes closed mitosis, where the nuclear envelope remains intact during cell division. Studying nucleoplasm properties over the course of closed mitosis may help reveal how nuclear volume, shape, surface area expansion, and chromosome segregation are linked and coordinated. To measure nucleoplasm material properties in S. pombe, we have paired Fluorescence Correlation Spectroscopy (FCS) with a photoswitchable fluorophore, enabling fine control over fluorescent intensity inside live cells. We infer material properties from FCS measurements, while the photoswitchable probe enables confocal imaging in conjunction with these measurements, yielding corresponding information about cellular state and dynamics. Interestingly, we find that nucleoplasm material properties do not vary significantly over the cell cycle. Future studies will use this tool to examine how diverse molecular and genetic perturbations alter nucleoplasmic properties, providing insight into how these properties maintain nuclear function and protect genomic integrity over the cell cycle and during development.

Restricting the localization of centromere-specific histone H3 variant Cse4 (CENP-A in humans) to centromeric chromatin is essential for chromosome segregation. Mislocalization of overexpressed Cse4/CENP-A to non-centromeric regions contributes to chromosomal instability (CIN) in model organisms and human cells. CIN is an important hallmark of many cancers and hence defining mechanisms that prevent mislocalization of Cse4 is clinically significant. Here we report a role for YRA2 (Yeast RNA Annealing Protein 2) in ubiquitin mediated proteolysis of Cse4 to prevent its mislocalization for chromosomal stability. YRA2 was identified in a genome-wide screen for gene deletions that exhibit synthetic dosage lethality (SDL) upon overexpression of CSE4 (GALCSE4). We determined that yra2{Delta} strains exhibit increased Cse4 stability, enriched Cse4 chromatin association, reduced Cse4 ubiquitination, Cse4 mislocalization, and CIN. Defects in interaction of E3 ubiquitin ligase Psh1 with Cse4 contributes to stability of Cse4 in yra2{Delta} strains. Consistent with these results, overexpression of PSH1 suppresses GALCSE4 SDL in yra2{Delta} strain. We determined that Yra2 mediated proteolysis of Cse4 is independent of its RNA related functions as strain deleted for the C-terminal ChTOP domain of Yra2 with an intact N-terminal RNA binding domain exhibits GALCSE4 SDL and defects in Cse4 proteolysis. Furthermore, poly(A)+ RNA export mutants in YRA1 (yra1-2) and MEX67 (mex67-5), that interact with Yra2, do not exhibit GALCSE4 SDL and defects in RNA export are not observed in yra2{Delta} cells. In summary, we have defined a key role for Yra2 in preventing mislocalization of Cse4 by facilitating its proteolysis to preserve chromosomal stability. Article summaryAccurate segregation of chromosomes during cell division is essential because segregation errors are linked to cancer and developmental disorders. We investigated how cells prevent mislocalization of centromere-specific histone H3 variant Cse4, which is essential for faithful chromosome segregation. We found that the yeast RNA annealing protein Yra2 prevents Cse4 mislocalization by promoting Psh1 mediated ubiquitination and degradation of Cse4. Cells lacking Yra2 showed increased stability of Cse4, enhanced chromatin enrichment with mislocalization to non-centromeric regions and CIN. These defects were suppressed by induction of Psh1. Our findings reveal a novel role for Yra2 in regulating Cse4 levels for chromosomal stability.

We present a chromosome-level genome assembly of a female Norwegian white-throated dipper (Cinclus cinclus cinclus) generated using Oxford Nanopore Technologies (ONT) long reads and Hi-C scaffolding. The assembly comprises two pseudo-haplotypes, hap1 (1186 Mb) and hap2 (1115 Mb), with 96.7% and 94.4% of sequences assigned to chromosome-scale scaffolds, respectively. Both pseudo-haplotypes contain 40 autosomes, with the Z and W sex chromosomes assigned to hap1. Compared with the PacBio HiFi-based C. c. gularis reference assembly bCinCin1.1.pri, which contains 38 autosomes, sequence represented as a single dot-chromosome (chr 36) is resolved into three distinct dot-chromosomes (chr 36, 39, and 40), a configuration supported by Hi-C contact patterns. BUSCO completeness was high for hap1 (99.2%) and hap2 (95.0%), with 19,003 and 17,746 predicted protein-coding genes, respectively. Compared with the HiFi-based C. c. gularis reference and HiFi-based assemblies generated from the same individual, the ONT-derived assemblies were substantially less fragmented and recovered more sequence from the smallest chromosomes. Synteny was otherwise largely conserved between subspecies. HiFi depletion increased strongly from macrochromosomes to micro- and dot-chromosomes, and HiFi-depleted regions were enriched for repeats and predicted non-B-DNA-associated features, particularly G-quadruplexes and direct repeats, whereas ONT coverage remained comparatively stable. These results show that conventional genome-wide assembly metrics can obscure substantial differences in the recovery of repeat-rich avian dot-chromosomes and highlight the value of chromosome-aware evaluation and ONT sequencing for recovering these regions.

Determining the correct donor species/lineages of grass-to-grass lateral gene transfer (LGT) is vital for deducing specific donor features that could help inform the mechanism of transfer. This requires a dataset spanning a broad range of species to achieve the phylogenetic resolution necessary for precise donor inference. As grass-to-grass LGT often involves the transfer of multi-gene DNA fragments, they can contain additional sequences that allow for accurate orthologous comparisons, such as nuclear DNA of plastid origin (NUPTs). Here we systematically scan for NUPTs in the genomes of four Alloteropsis semialata accessions, whose LGTs have previously been characterised. Using the abundant Panicoideae chloroplast sequences, we reconstruct NUPT phylogenies and infer two lateral acquisitions: one from Paniceae/Digitaria and another from Andropogoneae/Eremochloa adjacent to a previously identified LGT. We then assembled and included an additional 12 Eremochloa chloroplast genomes in the analysis and showed the likely donor was Eremochloa attenuata. Subsequent short-read mapping from E. attenuata to the nuclear region flanking this NUPT showed consistent coverage across the region, including the previously identified LGT, supporting co-transfer. Overall this study highlights the potential for NUPTs to better identify the donors of grass-to-grass LGT.

Telomeres are nucleoprotein structures that protect chromosome ends and are maintained by the Telomerase Reverse Transcriptase (TERT) protein that uses a noncoding Telomerase RNA (TR) as a template. In monkeyflowers, Mimulus lewisii had an ancient TR gene duplication, synthesizing an evolutionarily atypical sequence heterogeneous telomere. How TERT interacts with both TR paralogs during telomere maintenance is unknown and answers can shed novel insights underlying telomere function. Using new genome assemblies we discovered TERT is rapidly evolving in lineages sharing the TR duplication. We investigated the functional consequences arising from the rapid evolution, first by using yeast three-hybrid and testing the physical binding between conspecific and heterospecific TERT-TR combinations. Results showed TERT binds both ancestral (TR1) and derived (TR2) TR paralogs in M. lewisii, but not in species without a functioning TR2. We located the region of TR binding to amino acids near the KRxR motif. We then combined next-generation sequencing with Telomeric Repeat Amplification Protocol and discovered M. lewisii had high telomerase activity. Comparative transcriptomics indicated no strong evidence of expression divergence in telomere maintenance genes for M. lewisii, suggesting rapid evolution shaped TERT protein sequence. In vivo activity of M. lewisii telomerase was investigated by analyzing F1 telomeres generated by crossing M. lewisii and M. verbenaceus, which doesnt have a functioning TR2. Results showed M. verbenaceus chromosome ends in the F1 had converted into M. lewisii telomeres, suggesting dominance of the M. lewisii telomerase. We demonstrate TERT-TR coevolution can have significant consequences on the evolution of plant telomeres. Significance statementTelomeres protect chromosome ends and are maintained by the telomerase complex. We discovered the catalytic component of the telomerase (TERT) was rapidly evolving in monkeyflowers (Mimulus) and studied the molecular consequences. In M. lewisii, TERT evolved lineage-specific amino acids to bind two sequence divergent telomerase RNA paralogs. Telomerase activity assay showed M. lewisii synthesized more telomere repeats compared to its sister species without the TR duplication, and transcriptomics indicated this was not due to a change in telomere maintenance gene expression. Genetic experiments in interspecies hybrids showed M. lewisii telomerase could convert chromosome ends in sister species into M. lewisii-like telomeres suggesting functional dominance. We show rapid evolution of the telomerase can have significant effects on telomere evolution.

Complex rearrangements are one of the rarest types of structural variants (SVs) and can be divided into two categories: complex chromosomal rearrangements (CCRs) and complex genomic rearrangements (CGRs). CCRs include structural rearrangements that present at least three breakpoints and show exchange of genetic material between more than two chromosomes and CGRs are rearrangements that present more than one junction and/or more than one SV in cis. They are usually formed by one of the chromoanagenesis mechanisms, where a massive disruptive cellular event leads to multiple structural rearrangements. Classical cytogenomic techniques have been commonly applied for their characterization, but methodologies that involve longer DNA molecules, namely optical genome mapping (OGM) and long-read genome sequencing (lrGS), present a considerably higher SV detection resolution, revealing more details about the rearrangements, including precise breakpoint location. Here, we describe six patients with complex rearrangements investigated through a combination of different techniques: karyotyping, chromosomal microarray, and OGM were performed to characterize the rearrangements. Subsequently, lrGS was used to further resolve the alterations, refine their breakpoints' location, and sequence their junction points. Three patients presented CCRs involving three, four, and six chromosomes, while three exhibited CGRs involving one different chromosome each, providing a variety of complex SVs to show the importance of each technique and their combination in rearrangement resolution. In total, the complex rearrangements presented 127 breakpoints, 66 junction points and involved 14 of the 24 chromosomes. Higher-resolution techniques revealed additional complexity in all cases. Despite the advances provided by OGM and lrGS, conventional karyotyping remained indispensable for complete rearrangement resolution. In two patients, the findings supported a novel mechanism combining features of the different chromoanagenesis processes. Furthermore, evidence of inherited alterations was identified, and the comprehensive characterization of the rearrangements enabled more accurate genotype-phenotype correlations. Our findings indicate that an integrated approach combining karyotyping, OGM, and lrGS can completely resolve SVs, including complex rearrangements.

This paper presents and critically reviews the results of original and some literature based experimental studies conducted by the authors last years on the structural organization of DNA in dormant (starvation stress), anabiotic dormant (4 HR treatment) E. coli cells, as well as the K12 {Delta}dps strain, which lacks the Dps protein (Dps null E. coli). The experimental data includes small-angle synchrotron radiation diffraction (SAXS) and transmission electron microscopy (TEM) data. Synchrotron radiation diffraction experiments on K12{Delta}dps cells allowed us to conclude that peaks at 44.3, 22.1, and 14.8 angstrom resolutions are associated exclusively with ordered DNA organization. Peaks at 44.3, 22.1, and 14.8 angstrom resolutions are also observed for samples of dormant (starvation stress) cells and anabiotically dormant cells. Therefore, this ordered DNA organization also applies to samples of dormant and anabiotically dormant cells. A model is proposed that considers the ordered DNA organization in the cell as a cholesteric liquid crystal. The powder diffraction pattern calculated based on this model is compared with experimental small angle X ray scattering (SAXS) data obtained on Dps-null cell samples. The model completely reproduces the key features of the experimental diffraction pattern from Dps-null cell samples. Accordingly, the cholesteric liquid crystal model corresponds to DNA packaging in dormant and anabiotically dormant cells. Cholesteric liquid crystal ordering should be further considered in all models of cellular DNA packaging. To address the question of which structural organization of DNA predominates in the cell: the cholesteric liquid crystal or nanocrystalline or whether they coexist and fully manifest themselves under different external conditions, it is necessary to utilize the latest methodological advances in structural analysis.

Chromatin accessibility has long been used as a marker for regions of DNA with regulatory potential. Fiber-seq detects chromatin accessibility on individual DNA fibers, enabling analyses beyond the identification of the accessible chromatin regions (ACRs). By providing single molecule level high resolution, Fiber-seq provides unprecedented qualitative descriptions, including potential categorizations of ACRs, identification of internal transcription factor footprints and nucleosome positioning within individual DNA fibers. As with all tools, the power of this technique depends on careful experimental design and data analysis -- incorrect usage will result in incorrect conclusions. Here we offer guidelines and flag potential pitfalls when generating and analyzing Fiber-seq data, such as (1) the optimum levels of adenosine methylation per-fiber, (2) the power of per-fiber state inference, (3) the importance of controlling for read depth and methylation rates when comparing across samples, (4) the limitations of long-read sequence mapping, and (5) suggestions for identification of differentially accessible peaks across samples.

HP1 proteins are highly evolutionarily conserved chromatin-associated factors known to play essential roles in genome stability and nuclear organization. In mammals, three HP1 isoforms, HP1, HP1{beta} and HP1{gamma}, have been described, but their individual functions remain incompletely characterized. Here, we inactivated HP1 or HP1{beta} in different cell lines and quantified chromosomal breaks on metaphase spreads in the presence or absence of aphidicolin-induced replication stress. Loss of HP1, but not of HP1{beta}, led to a significant increase of chromosomal breaks on chromosome arms and within pericentromeric heterochromatin under these conditions. Mechanistically, loss of HP1 was associated with a reduction in replication fork velocity, suggesting that HP1 deficiency induces a replication stress that sensitizes specific genomic loci to replication perturbation. Consistent with this, HP1 loss was associated with a moderate but consistent increase in {gamma}H2AX and 53BP1 foci, an increased occurrence of DNA synthesis during mitosis, and enhanced recruitment of FANCD2, all recognized as hallmarks of common fragile site (CFS) expression. In addition, rescue experiments using a chromodomain mutant HP1 (V22M) unable to bind H3K9me3 indicated that HP1 protective function over these specific foci did not require its interaction with this histone mark. Altogether, these data indicate that, independently of its binding to H3K9me3, HP1 stabilizes specific genomic regions that behave as HP1-dependent fragile sites, at least in part by regulating replication fork progression, limiting mitotic DNA synthesis possibly by competing with FANCD2 for chromatin access at these regions.

The European leaf-toed gecko (Euleptes europaea) is a small, nocturnal gecko endemic to the western Mediterranean. As a phylogenetically distinctive member of the Gondwanan family Sphaerodactylidae, it represents an important species for studying Mediterranean island biogeography, adaptation, and reptile genome evolution. The species also occupies a key position for investigating the evolution of sex chromosomes, as geckos exhibit remarkable diversity and frequent transitions in sex-determination systems. We present a chromosome-level genome assembly of Euleptes europaea generated as part of the Vertebrate Genomes Project. The 1.8 Gb assembly has a scaffold N50 of 102.3 Mb (contig N50 27 Mb), with 21 chromosome-scale scaffolds corresponding to the known karyotype (2n = 42). The primary assembly has a BUSCO completeness of 97.80% (95.60% as single-copy), a k-mer completeness of 96.00%, and a k-mer quality value (QV) of 61.20. Repetitive elements account for 53.20% of the genome and genome annotation identified 18,633 protein-coding genes. This high-quality reference genome will facilitate studies of genome evolution, island adaptation, and sex chromosome evolution across geckos and other reptiles.

Heterochromatin replication isnt random or uniform, but occurs in a characteristic spatial and temporal pattern. Previous studies produced conflicting models for pericentric heterochromatin (PCH) replication, suggesting either that heterochromatic sequences translocate to the domain surface for replication, or that replication can also occur internally through localized decondensation. To distinguish true overlap from peripheral enrichment around an irregular PCH domain, we developed FOC-Map, a colocalization analysis approach that combines segmentation of one channel with binning of the other. Applying FOC-Map to three-dimensional Airyscan imaging of cultured Drosophila cells, we find that replication foci at the onset of late S-phase are confined to the outer boundary of the PCH domain, forming a shell-like pattern with little overlap into the HP1a-rich interior. As late S-phase progresses, replication foci are observed within the domain, localizing to low-HP1a regions interspersed between more condensed regions. We then assessed the distribution of CDC45, a rate-limiting replication initiation factor, and found that CDC45 foci are depleted from the PCH domain throughout the cell cycle. We propose that low levels of CDC45 within HP1a-rich PCH limit replication initiation to the domain periphery, giving rise to the shell-like pattern of replication foci that progressively works inward until PCH replication is complete.

BackgroundLamins are major regulators of the spatial and functional organization of chromatin. Lamins at the nuclear periphery form the lamina that anchors heterochromatin to the nuclear envelope. A subpool of A-type lamins localizes in the nuclear interior, where they also bind to euchromatic genomic regions. A-type lamin properties and chromatin association are regulated by lamin-associated polypeptide 2alpha (LAP2). Here we systematically analyze, how LAP2 depletion affects chromatin organization, accessibility and gene expression on a genome-wide level. ResultsLAP2 depletion in mouse dermal fibroblasts positively and negatively affects chromatin accessibility and gene expression throughout the genome, which correlates with changes in chromatin association of A-type lamins and the nucleosomal remodeler proteins BRG1 and CHD4. In particular, A-type lamins bind to open chromatin regions close to BRG1 and CHD4 binding sites and deregulated genes, but do not directly accumulate on genes and BRG1 and CHD4-enriched sites. Unsupervised clustering of the datasets on LAP2-bound genomic regions confirms spreading of A-type lamins to active chromatin regions containing deregulated genes and an enrichment of chromatin remodelers on a subset of these genomic regions. ConclusionsLAP2 depletion in fibroblasts leads to a gross rearrangement of chromatin. Genome-wide chromatin reorganization is linked to spreading of A-type lamins to active chromatin regions and accompanied by a restriction of chromatin remodelers to a subset of active genomic regions. These changes correlate with changes in chromatin accessibility and gene expression throughout the genome, particularly in regions where lamin binding is gained in LAP2 knockout versus wildtype cells.

Damage to biological specimens by the electron beam is the fundamental resolution-limiting factor in cryoelectron microscopy (cryo-EM) single particle analysis. There is, however, currently no method to accurately infer fluence-dependent changes to the specimen structure during electron irradiation. We develop a Bayesian framework to fit a sequence of atomic models to a series of cryo-EM reconstructions produced at increasing fluence. In particular, our algorithm is able to infer the ensemble average position and atomic displacement parameter of every atom in the macromolecule as a function of fluence. Application of the algorithm to cryo-EM datasets shows that the molecule expands during imaging and identifies environment-dependent variations in beam-induced damage. We use our results to propose a stochastic process model of this phenomenon. We envisage that our method will lead to a better mechanistic understanding of radiation damage to biological specimens and may contribute to efforts to mitigate its effects.